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Abstract 

Time evolution of an advection-dominated accretion flow is explored in terms of relativistic fluid equa- 
tions. An axially symmetric, vertically averaged mean flow is constructed and then perturbed. An axisym- 
metric pattern is followed as it propagates in the form of a wave towards the horizon of a rotating (Kerr) 
black hole. Several assumptions are relaxed in comparison with previous works (Manmoto et al. 1996). 
The wave reflection and steepening are examined and the influence of black hole angular momentum is 
discussed. 
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1. Introduction 

In this paper we report on computations in which rel- 
ativistic hydrodynamics has been employed to describe 
propagation of a wave pattern in accreted fluid. An 
axisymmetric perturbation is imposed on a background 
model of a transonic, advection-dominated flow. The fluid 
has negligible self-gravity; it moves in spacetime of a Kerr 
black hole that is described by parameters M and a.0 
While our approach is essentially one-dimensional, it is 
able to treat non-linear regime of the growing perturba- 
tion. The model can provide some insight into phenomena 
expected to occur in astrophysically realistic models. 

We build our computations on the results from several 
earlier papers. Particularly relevant have been computa- 
tions of Manmoto et al. (1996; cited as PI hereafter) who 
explored wave propagation and reflection in a transonic 
flow (see also Kato et al. 1998, chapt. 11). These authors 
found clear evidence for reflecting waves near the hori- 
zon. The perturbation then steepens into a shock wave. 
Observational motivation for these studies comes from the 
evidence for X-ray fluctuations, origin of which is often 
linked with disturbances in the innermost regions of ac- 
cretion flows. In PI, the background accretion flow was 
constructed in the pseudo-Newtonian framework, while 
we employ general relativistic description and relax sev- 
eral other restrictions. 

The main aim of this note is to show the results of our 
computations, where black hole rotation stands as one of 
the parameters. We refer to other, more extended pa- 
pers for a thorough discussion of various connections and 
the place that these computations may have in a general 
scheme of the black-hole accretion. For instance, a lot of 



Hereafter, equations will be written in gcometrizcd units and 
scaled with the central mass M . Hence, dimensionless angular 
momentum acquires values 0<<i < 1- Corresponding quantities 
in physical units can be obtained by straightforward conversions 
(Misner et al. 1993, p. 36). Gravitational radius, ijg = GM j (? , 
will be adopted as a convenient unit of length. 



interest has been focused on the possibility of limit-cycle 
oscillations in accretion disks (Honma et al. 1991). Indeed, 
non-linear time-dependent calculations of vertically inte- 
grated transonic accretion disks exhibit limit-cycle be- 
haviour associated with thermal instability (Szuszkiewicz 
& Miller 1998). Recently, Bate et al. (2002) analyzed ax- 
isymmetric waves with a two-dimensional hydrodynamical 

scheme. 

This paper is organized as follows. In Sec. 2T the back- 
ground stationary flow is constructed. 
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In Sec. 

turbation is imposed on the background flow and its time 
evolution is examined. Finally, the properties of the tem- 
poral behaviour are summarized in Sec. |3| 

2. Equations and the method of solution 

We assume Kerr black-hole spacetime^ in which hydro- 
dynamic equations are solved. For the background flow, 
we adopt a model described by Gammie & Popham (1998) 
and Popham & Gammie (1998). Steady-state equations 
are summarized in subsection for convenience, while 
our main attention will be focused on time evolution of a 
perturbation discussed later on in the text. 

2.1. A steady background accretion flow 

The steady-state equations can be derived from: 

(i) Particle conservation {pu^).^ = which, thanks to 
stationarity, axial symmetry and vertical averaging 
can be expressed in terms of accretion rate: M — 

(ii) Energy conservation u^T'^'^-i, ~ 0, which can be writ- 
ten in the form u^u,r — {u + p)p~-^u^ p^r = $ — A=/$, 
where $ = —t^^(T^^\ f represents the ratio of heating 
and cooling processes, respectively. 



We use Boyer-Lindquist coordinates {t,r,6,<l); Misner et al. 1973) 
with usual notation; T, = + a'^ cos^ 9 , A = r^ — 2r + a^, A={r^ + 
a2)2- Aa^sin^e, = l + a^ /r^ + 2a'^ /r^ , = 1- 3/r + 2a/r'^^^ , 
& = 1 — 2/r + /r^; g = siv? 9 is determinant of the metric. 
Horizon is at r.)- = 1 -|- — d^, 9 = tx /2 is the flow mid-plane. 
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Fig. 1. Time evolution of surface density perturbation in our 
model (5S in units of g/cm'^). Schwarzschild metric was as- 
sumed with a = 0.1, / = 1- 

(iii) Radial component of the momentum conservation 
law, hr^Tt'^.^^O. 

(iv) Angular momentum conservation, T^^S^f^.^y = 0. 

In analogy with (Abramowicz et al. 1997), we defined the 
aspect ratio of the flow, Hq =H/r (where H is geometrical 
semi-thickness); h^p denotes the projection tensor, is 
the axial Killing vector, and a^i, is the equilibrium shear 
stress tensor. 

Four equations correspond to conservation laws (i)-(iv): 

(1) 
(2) 

(3) 
(4) 
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The unknowns are: mid-plane density and temperature, 
p and T, and two components of the fluid four-velocity 
In eq. (|), we denoted jT = 1 - 2ar2 - (r^ - a^) Vt^ , 
VI = v!^ jv}' . In eq. (Q), j is the eigenvalue of the problem 
and is the viscous tensor component. See Gammie 
& Popham (1998) and the Appendix for details of the 
derivation. 

Equations (|l|)-(^ were solved numerically and the ob- 
tained results agree with Gammie & Popham (1998). The 
steady solution will be used as the initial state of the time- 
dependent flow integration in the next subsection. 

2.2. Time evolution of a perturbed flow 

Equations describing evolution of the flow are equiva- 
lent to those listed in items (i)-(iv) in the previous sub- 
section, except that we do not neglect terms containing 
time derivatives. 



^ One of the two components can be chosen in the azimuthal di- 
rection, l = u^, the latter is conveniently expressed as the radial 
speed V of the fluid local rest frame (LRF) with respect to the 
local frame co-rotating with the fluid at constant radius (CRF). 
In eq. (p|), 7^ is the Lorentz factor corresponding to azimuthal 
boost between CRF and the locally non-rotating frame (LNRF) 
and similarly 7^ = (1 — V'^)^^/'^ . 



The initial imposed perturbation is adopted analogi- 
cal as in ref. PI (sec. 2), so that comparisons are eas- 
ier. Namely, for the density perturbation we assume 
5p/p — Kgcxp [— (r — ro)^/(Ao/2)^] , where kq < 1 is the 
starting perturbation amplitude, tq is its location and Aq 
is the characteristic length. Notice that density perturba- 
tion can be related to perturbations of other quantities - 
T, /, and V (Kato et al. 1996). An explicit form of the 
equations is given in the Appendix. We solved the set 
(Al)-(A5) using the Lax-WendrofF scheme (Press et al. 
1992) with the free ingoing boundary condition imposed 
close above r = r_|_ (below the sound point). Far from the 
black hole, the asymptotic form of the flow goes over to 
the stationary solution, eqs. (0)-(U). 

3. Results: Waves propagation and reflection 

First, we reproduce (in the general relativistic frame- 
work) the basic behaviour of the flow perturbation. 
Namely, we can observe the initial propagation of the wave 
towards horizon, its reflection in the region of several i?g 
(near the sonic point) and then formation of a shock. 

In Figure |l|, perturbation (5I](r, t) of the surface den- 
sity E 2rH0p is plotted. We set M = IQMq and 
M = 0.16Le/c2 = 2.24 x 10"g s"! for definiteness. The in- 
tegration domain was 300i?g x 10^i?g/c (somewhat smaller 
extent is captured in the graph; Rg/c corresponds to 
~ IQ-'^s for M = 10^0 ). This figure can be compared 
with Fig. 1 in PI (after a trivial scaling of units). Our 
integration time is longer, allowing us to observe subse- 
quent outgoing waves appearing at t > 3000 (« 10 viscous 
periods at the outer edge of the integration grid). For 
definiteness we adopted a = here. 

In Figure 1^, 5Yi{r,t) is shown for two limiting values 
of angular momentum, a ~ Q and a = 1, respectively. 
Contour-lines show a crest of the initial perturbation 
which proceeds with the accretion flow. At small radii, 
pattern speed is somewhat influenced by black hole an- 
gular momentum. Several reflected waves appear at later 
stages until the simulation settles back to the initial state. 
Given the initial perturbation and keeping a and / fixed, 
the form and the amplitude of waves are slightly sensitive 
to the black hole angular momentum - the effect that is 
difficult to treat within the pseudo-Newtonian approach. 
However, the general scheme is in good agreement with 
the pseudo-Newtonian treatment of paper PI. 

4. Conclusions 

We examined temporal behaviour of a simple relativis- 
tic model of the advection-dominated flow. To this aim, 
we observed the evolution of disturbances originating in 
outer parts of the disk. Reflection of the waves was dis- 
cussed previously in various contexts. We developed a 
robust code which is able to deal with the problem across 
the large range of disk parameters without introducing 
numerical viscosity. The scheme appears useful as a test 
bed for more complicated investigations with unstable be- 
haviour, treated in more dimensions (work in progress). 
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Fig. 2. Surface density perturbation for a = (left) and a = 1 (right). The perturbation is imposed at r = 100 (time t = 0) and it 
proceeds towards the black hole. The first outgoing wave is reflected at about t = 1500 where contour-lines indicate the maximum 
perturbation. Further waves appear at later times. Contours are plotted with equidistant spacing. Parameters a and / as in Fig. hi. 



In this note we described a stable system which after 
a transition period eventually relaxes to the initial steady 
state. Our scheme, thanks to its rather general formula- 
tion, exhibits more complex behaviour of reflected waves 
than the pseudo-Newtonian model described in PI. The 
main advantages are general relativistic treatment of hy- 
drodynamic equations including finite propagation speed 
of viscous effects. We are also able to follow the evolu- 
tion for longer intervals, till the shock escapes from the 
observed region and the model settles back to its initial 
equilibrium state. 

VK acknowledges helpful discussion with Professor 
Marck Abramowicz. Support from the grants GACR 
205/00/1685 and 202/02/0735 is also acknowledged. 

Appendix. 

The only non- vanishing component of the viscous tensor 
(in LRF), S = 't{r){ct>)^ satisfies equation u^^S-^ = —t~^(S — 
Sq) where is the fluid relaxation time. An equilib- 
rium value Sq is given in terms of a parameter, speed of 
sound Cg, enthalpy rj = {p-\-u+p)/ p, and the shear tensor 
component a = a^r){4>)- Sq = —2pr]i/a, where = aCsrHg. 
Energy conservation can be further adapted using equa- 
tion of state, p = pT, where p is the rest mass density and 
T is dimensionless temperature (T = kT]^/rnc^ where Tk is 
kinetic temperature and to is the mean molecular mass). 



We summarize time-dependent equations of sec. |2.2 
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Here, the unknown functions (T, I, V, p, S) all depend 
on t and r (variables I and V were introduced in foot- 
note 3). Geometrical semi-thickness of the flow is given 
by p7] — {u\ — l)] Hi = pr^. Time dependence of Hg 
can be expressed via T, I and Y . Finally, for the dissipa- 
tion function one flnds $ = —2Sa(r){^) with <J(^r){4>) being 
the LRF time-dependent shear component, which is of the 
same form as in steady state (do) plus correction terms: 



^2 ' dt 2rV^df 



(A6) 
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